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ABSTRACT 

We  present  orbital  parameters  for  six  double-lined  spectroscopic  binaries  (i  Pegasi,  co  Draconis,  12  Bootis,  VI 143 
Cygni,  ft  Aurigae,  and  Mizar  A)  and  two  double-lined  triple  star  systems  ( k  Pegasi  and  rj  Virginis).  The  orbital 
fits  are  based  upon  high-precision  radial  velocity  (RV)  observations  made  with  a  dispersed  Fourier  Transform 
Spectrograph,  or  dFTS,  a  new  instrument  that  combines  interferometric  and  dispersive  elements.  For  some  of  the 
double-lined  binaries  with  known  inclination  angles,  the  quality  of  our  RV  data  permits  us  to  determine  the  masses 
Mi  and  M2  of  the  stellar  components  with  relative  errors  as  small  as  0.2%. 

Key  words:  binaries:  spectroscopic  -  instrumentation:  spectrographs  -  techniques:  radial  velocities 


1.  INTRODUCTION 

For  the  past  several  years,  our  research  group  has  been  devel¬ 
oping  a  new  optical  spectrograph  concept  called  the  dispersed 
Fourier  Transform  Spectrograph,  or  dFTS.  The  instrument  de¬ 
sign  merges  a  traditional  Fourier  Transform  Spectrometer  with 
a  dispersive  grating  spectrograph  such  that  the  interferometer 
output  is  divided  into  thousands  of  narrowband  channels,  all 
operating  in  parallel.  This  multiplex  advantage  boosts  the  effec¬ 
tive  throughput  of  the  system  by  a  large  factor,  making  the  dFTS 
competitive  with  echelle  spectrographs  for  spectroscopic  anal¬ 
ysis  of  stars,  particularly  measurement  of  their  radial  velocities 
(RVs). 

Hajian  et  al.  (2007)  describe  our  prototype  device,  dFTSl, 
and  explain  the  underlying  theory  and  hardware  implementation 
in  detail.  Based  upon  our  commissioning  observations  with 
dFTS  1 ,  we  subsequently  designed  and  built  a  second-generation 
version,  dFTS2,  which  we  deployed  to  the  Steward  Observatory 
2.3  m  Bok  Telescope  for  a  year-long  observing  campaign.  In 
Behr  et  al.  (2009),  we  discuss  the  dFTS2  hardware  and  present 
velocimetry  measurements  of  RV  standard  stars  and  single-lined 
spectroscopic  binary  stars  (SBls). 

In  this  paper,  we  describe  the  results  from  our  dFTS2 
observations  of  double-lined  spectroscopic  binaries  (SB2s)  and 
double-lined  triple  systems.  SB2s  provide  one  of  the  best  means 
for  measuring  the  masses  of  stars:  given  an  accurate  RV  curve 
for  each  stellar  component  and  the  inclination  angle  i  of  the 
orbital  plane  to  the  observer’s  line  of  sight,  we  can  derive 
the  component  masses  using  Kepler’s  third  law.  Traditional 
spectroscopic  observations  with  an  echelle  spectrograph  and 
a  thorium-argon  calibration  source  can  achieve  a  velocity 
precision  of  ~0.01-0. 10  km  s~ 1  on  late-type  narrow-lined  stars 
(Ramm  et  al.  2004;  Skuljan  et  al.  2004;  Tomkin  &  Fekel 
2006;  Fekel  et  al.  2007;  Ramm  2008).  For  greater  precision, 
Konacki  (2005, 2009)  has  developed  a  technique  using  an  iodine 
absorption  cell  with  which  the  RVs  of  a  spectroscopic  binary 
can  be  measured  at  the  0.005-0.010  km  s_1  level.  Our  dFTS2 
instrument,  in  contrast,  achieves  high  RV  precision  and  stability 
without  a  superposed  reference  spectrum.  As  described  in  Behr 


et  al.  (2009),  we  measure  the  RVs  of  non-binary  stars  and  single- 
lined  binaries  to  0.01-0.03  km  s~  1  and  anticipate  even  better 
performance  once  thermal  stability  issues  in  our  instrument 
design  have  been  addressed. 

2.  DATA  ACQUISITION  AND  RV  ANALYSIS  PROCEDURE 

The  data  reported  in  this  paper  were  collected  between 
2007  October  and  2008  June  during  bright-time  observing  runs 
at  the  2.3  m  Bok  Telescope  of  the  Steward  Observatory  on 
Kitt  Peak.  An  observation  of  a  given  SB2  target  consisted 
of  500  exposures  spanning  a  range  of  interferometer  delays 
corresponding  to  a  spectral  resolution  of  approximately  50,000. 
Each  exposure  lasted  1. 0-4.0  s  in  duration,  depending  on  the 
star’s  brightness  and  the  atmospheric  seeing  and  opacity.  Each 
scan  also  required  a  total  overhead  time  of  approximately  four 
minutes,  independent  of  exposure  time,  for  CCD  readout  and 
moving  to  the  next  delay  position.  Our  targets,  listed  in  Table  1, 
were  chosen  because  of  their  relative  brightness  and  short 
periods,  so  that  we  could  acquire  many  observations  per  star 
during  the  limited  duration  of  this  initial  observing  campaign. 

To  measure  double-lined  RVs  from  our  interferogram  data, 
we  employed  a  variation  of  the  standard  two-dimensional 
cross-correlation  technique  (Mazeh  &  Zucker  1994).  Instead 
of  transforming  our  interferograms  into  spectra  and  then  cross- 
correlating  model  templates  against  each  observed  spectrum, 
we  converted  the  template  spectra  into  template  interferograms 
and  then  compared  those  model  interferograms  to  the  observed 
interferograms.  Because  interferograms  add  linearly,  we  can 
compute  sequences  of  single-lined  interferograms  for  the  A 
and  B  components  of  a  binary,  spanning  a  range  of  RVs  for 
each,  and  then  add  together  A  and  B  interferograms  to  create 
a  two-dimensional  grid  of  double-lined  model  interferograms. 
Calculating  the  x2  difference  between  the  models  and  the 
observed  interferogram  data,  we  construct  a  map  of  fit  quality, 
where  the  minimum  point  indicates  the  best-fit  solution  for 
component  velocities  Vi  and  V  2,  and  the  projection  of  the 
X2  =  Xmin  +  2.30  contour  line  onto  the  V 1  and  V 2  axes  provides 
la  error  bars  on  each  RV  point.  The  shapes  of  the  x2  contours 
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Table  1 

Spectroscopic  Binary  Targets  Observed  with  dFTS2 


Star 

V  Magnitude 

Spectral  Type 

P 

(days) 

V  sin  iA,B 
(km  s-1) 

References 

(  Peg 

3.8 

F5V 

10.21 

7.6, 7.2 

Fekel  &  Tomkin  (1983) 

co  Dra 

4.8 

F5V 

5.28 

7.1,  6.6 

Mayor  &  Mazeh  (1987),  Fekel  et  al.  (2009) 

12  Boo 

4.8 

F8IV 

9.60 

13.1,  10.4 

Boden  et  al.  (2005),  Tomkin  &  Fekel  (2006) 

VI 143  Cyg 

5.9 

F5V 

7.64 

23.6,36.2 

Andersen  et  al.  (1987) 

fS  Aur 

1.9 

A2IV 

3.96 

34.5,  35.0 

Smith  (1948),  Pourbaix  (2000) 

Mizar  A 

2.3 

A2V 

20.54 

32.6,36.2 

Fehrenbach  &  Prevot  ( 1961),  Pourbaix  (2000) 

k  Peg 

4.2 

F5IV 

5.97 

7.1,47.9 

Mayor  &  Mazeh  (1987),  Flajian  et  al.  (2007) 

r]  Vir 

3.9 

A2V 

71.79 

5.1,  4.5 

Hartkopf  et  al.  (1992),  Hummel  et  al.  (2003) 

were  elliptical  with  minor  and  major  axes  aligned  to  the  V\ 
and  V 2  axes,  indicating  no  significant  covariance  between  the 
component  velocity  errors. 

For  optimal  results,  the  template  spectra  must  be  well 
matched  to  the  actual  spectra  of  the  two  stellar  components. 
We  generated  synthetic  spectra  using  the  SPECTRUM  spec¬ 
tral  synthesis  package  (Gray  &  Corbally  1994;  see  also  http:// 
www.phys.appstate.edu/spectrum/spectrum.html)  and  then  var¬ 
ied  the  transition  strength  log  gf  of  each  line  and  the  projected 
rotation  velocity  v  sin ;  for  each  stellar  component  to  minimize 
the  x2  difference  between  the  model  and  the  observed  data.  A  fi¬ 
nal  template  spectrum  for  each  component  was  calculated  from 
a  median-filtered  average  of  transition  strengths  from  the  in¬ 
dividual  observations,  and  the  final  two-dimensional  RV  cross¬ 
correlation  was  then  performed  using  these  templates.  We  found 
that  we  could  derive  relatively  precise  and  self-consistent  v  sin  i 
values  from  our  interferograms  (as  listed  in  Table  1);  a  future  pa¬ 
per  will  explore  the  use  of  dFTS  data  to  measure  stellar  rotation 
velocities  and  other  line  broadening  mechanisms. 

Because  our  template  spectra  are  generated  using  atomic  tran¬ 
sition  wavelengths  from  the  National  Institute  of  Standards  and 
Technology  catalog,  the  derived  RVs  can  be  considered  “accu¬ 
rate”  in  the  sense  that  they  reflect  the  total  Doppler  shift  between 
the  rest  wavelength  of  a  line  and  the  observed  wavelength.  We 
have  not  (yet)  attempted  to  tie  our  velocity  scale  to  any  IAU 
velocity  standards,  nor  do  we  make  any  correction  for  gravi¬ 
tational  redshift  effects.  The  only  adjustment  made  to  the  RV 
data  is  conversion  to  a  solar  system  barycenter  reference  frame, 
using  the  IRAF  tool  bcvcorr.  These  barycentric  RV  data  are 
listed  in  Table  2.  It  should  be  noted  that  for  the  k  Peg  RVs,  “V i” 
refers  to  the  Bb  component  and  “  V2”  refers  to  the  A  component. 

We  derived  orbital  parameters  from  our  RV  data  points  using 
the  IDL  routines  CURVEFIT,  a  gradient-expansion  nonlinear 
least-squares  fitting  algorithm  included  with  the  IDL  package, 
and  HELI0_RV  (Landsman  1993),  which  computes  a  line- 
of-sight  velocity  curve  for  a  binary  component  given  the 
period  P,  periastron  time  T  (or  for  circular  orbits,  the  time  of 
maximum  positive  velocity),  eccentricity  e,  periastron  longitude 
co,  RV  semi-amplitude  K,  and  systemic  velocity  Vo  (alternatively 
denoted  as  y  by  some  researchers).  We  fit  the  primary  and 
secondary  RV  points  simultaneously,  assuming  that  P,  T,  e, 
and  Vo  are  the  same  for  both  components  and  that  co\  and  ah 
differ  by  180°.  The  CURVEFIT  routine  returns  1  a  uncertainties 
(standard  deviations)  for  all  derived  parameters.  For  all  of  our 
SB2  targets,  we  adopted  the  orbital  period  P  from  previously 
published  analyses  because  our  observations  covered  a  relatively 
short  period  of  time.  We  did  not  correct  for  the  light  travel  time 
across  each  binary  system,  because  the  resulting  changes  in 
the  RV  values  are  small  compared  to  the  RV  error  bars  in  all 
six  cases. 


Figure  1.  RV  measurements  of  the  SB2  system  1  Pegasi.  Filled  squares  show 
the  measured  RV  of  the  primary  component,  and  open  squares  indicate  the 
secondary  component.  Our  observing  campaign  had  ended  before  we  were  able 
to  complete  the  phase  coverage  of  this  system. 


3.  RESULTS  FOR  DOUBLE-LINED  BINARY  SYSTEMS 
3.1.  l  Pegasi 

Our  RV  measurements  for  the  double-lined  spectroscopic 
binary  1  Pegasi  (HR  8430,  HD  210027,  and  HIP  109176)  are 
plotted  in  Figure  1.  The  most  recent  published  RV  work  on 
this  system  comes  from  Fekel  &  Tomkin  (1983)  whose  orbital 
parameters  are  listed  in  Table  3  along  with  the  values  that  we 
derive  from  our  dFTS2  observations.  In  addition  to  adopting 
their  value  of  the  system’s  orbital  period,  we  also  followed 
their  lead  in  assuming  a  circular  orbit,  because  our  RV  points 
only  covered  half  of  the  orbital  phase,  and  CURVEFIT  could 
not  place  meaningful  constraints  on  e  or  co.  Our  values  for  K\ 
and  K2  are  compatible  with  those  of  Fekel  &  Tomkin,  although 
our  solution  for  Vo  differs  by  a  statistically  significant  amount. 
This  discrepancy  may  indicate  the  gravitational  influence  of 
an  unseen  and  distant  third  stellar  component  of  the  system, 
although  Tokovinin  et  al.  (2006)  did  not  find  any  close  tertiary 
companions  in  Two  Micron  All  Sky  Survey  images  of  1  Peg, 
and  the  astrometric  observations  of  Boden  et  al.  (1999)  saw  no 
evidence  for  a  companion  either.  Alternatively,  the  difference  in 
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Table  2 

RV  Data  Measured  with  dFTS2 


Star 

HJD  -  2.400,000 

Phase 

Vj  (km  s  1 ) 

V\  Error  (km  s  1 ) 

V2  (km  s  *) 

V2  Error  (km  s  1 ) 

t  Peg 

54,400.6388 

0.109 

33.302 

0.017 

-64.348 

0.080 

t  Peg 

54,401.7122 

0.214 

6.584 

0.019 

-21.902 

0.091 

i  Peg 

54,401.7269 

0.215 

6.235 

0.017 

-20.958 

0.084 

«  Peg 

54,402.6529 

0.306 

-20.809 

0.021 

22.781 

0.109 

«  Peg 

54,403.6604 

0.404 

-44.161 

0.020 

59.933 

0.103 

(  Peg 

54,404.6147 

0.498 

-52.709 

0.032 

73.279 

0.196 

t  Peg 

54,604.9660 

0.115 

31.993 

0.017 

-62.573 

0.080 

«  Peg 

54,606.9635 

0.311 

-22.230 

0.015 

24.831 

0.067 

«  Peg 

54,634.8758 

0.044 

42.325 

0.015 

-78.934 

0.067 

t  Peg 

54,634.8959 

0.046 

42.119 

0.016 

-78.687 

0.072 

t  Peg 

54,635.9132 

0.145 

25.341 

0.017 

-51.490 

0.076 

t  Peg 

54,635.9342 

0.147 

24.823 

0.016 

-50.598 

0.075 

t  Peg 

54,637.8866 

0.339 

-29.723 

0.017 

36.656 

0.073 

t  Peg 

54,637.9068 

0.341 

-30.242 

0.017 

37.434 

0.080 

t  Peg 

54,638.9426 

0.442 

-49.514 

0.015 

68.427 

0.066 

t  Peg 

54,638.9628 

0.444 

-49.713 

0.015 

68.752 

0.064 

w  Dra 

54,547.9814 

0.160 

-48.332 

0.027 

29.389 

0.049 

co  Dra 

54,549.9958 

0.542 

19.025 

0.026 

-53.678 

0.046 

co  Dra 

54,578.9388 

0.024 

-44.148 

0.031 

24.305 

0.056 

co  Dra 

54,579.9848 

0.222 

-41.842 

0.027 

21.495 

0.050 

co  Dra 

54,603.8020 

0.733 

13.203 

0.024 

-46.284 

0.045 

co  Dra 

54,603.9267 

0.756 

9.247 

0.029 

-41.483 

0.054 

co  Dra 

54,604.9389 

0.948 

-31.906 

0.026 

9.346 

0.049 

co  Dra 

54,605.8126 

0.114 

-49.790 

0.023 

31.110 

0.041 

co  Dra 

54,605.9378 

0.137 

-49.461 

0.026 

30.731 

0.048 

co  Dra 

54,634.9202 

0.627 

22.606 

0.025 

-57.902 

0.046 

co  Dra 

54,634.9758 

0.637 

22.298 

0.023 

-57.528 

0.040 

co  Dra 

54,635.6924 

0.773 

6.250 

0.022 

-37.884 

0.039 

co  Dra 

54,636.6824 

0.960 

-34.305 

0.024 

12.288 

0.042 

co  Dra 

54,636.9410 

0.009 

-42.260 

0.022 

22.042 

0.039 

co  Dra 

54,637.6927 

0.152 

-48.736 

0.027 

29.976 

0.049 

co  Dra 

54,637.9312 

0.197 

-45.017 

0.026 

25.330 

0.047 

12  Boo 

54,547.9153 

0.591 

-29.048 

0.065 

50.273 

0.072 

12  Boo 

54,548.9127 

0.694 

-48.375 

0.076 

70.324 

0.083 

12  Boo 

54,549.9111 

0.798 

-51.801 

0.080 

73.760 

0.093 

12  Boo 

54,577.7996 

0.702 

-49.395 

0.049 

71.239 

0.056 

12  Boo 

54,579.8185 

0.912 

-18.824 

0.060 

39.831 

0.070 

12  Boo 

54,579.8722 

0.918 

-16.030 

0.054 

37.061 

0.062 

12  Boo 

54,602.8623 

0.312 

50.008 

0.070 

-31.349 

0.079 

12  Boo 

54,603.7304 

0.402 

22.867 

0.073 

-3.201 

0.084 

12  Boo 

54,604.7012 

0.503 

-6.315 

0.061 

26.942 

0.071 

12  Boo 

54,605.7012 

0.607 

-32.713 

0.055 

53.903 

0.065 

12  Boo 

54,606.6579 

0.707 

-49.876 

0.063 

71.703 

0.071 

12  Boo 

54,635.7226 

0.733 

-52.290 

0.064 

74.099 

0.073 

12  Boo 

54,638.7263 

0.046 

58.158 

0.084 

-39.386 

0.092 

VI 143  Cyg 

54,603.8512 

0.742 

6.959 

0.260 

-39.988 

0.490 

VI 143  Cyg 

54,604.8585 

0.873 

57.691 

0.200 

-92.635 

0.333 

VI 143  Cyg 

54,605.8830 

0.007 

59.725 

0.191 

-95.517 

0.309 

VI 143  Cyg 

54,606.8736 

0.137 

-69.958 

0.250 

37.136 

0.445 

VI 143  Cyg 

54,634.8268 

0.796 

22.287 

0.168 

-55.975 

0.320 

VI 143  Cyg 

54,636.8467 

0.060 

-29.744 

0.205 

-0.671 

0.339 

VI 143  Cyg 

54,638.8573 

0.323 

-65.668 

0.192 

33.420 

0.335 

P  Aur 

54,400.9449 

0.134 

54.375 

0.151 

-91.605 

0.149 

P  Aur 

54,401.9169 

0.379 

-96.561 

0.158 

63.013 

0.156 

P  Aur 

54,402.9819 

0.648 

-81.719 

0.193 

48.833 

0.195 

P  Aur 

54,402.9909 

0.651 

-81.088 

0.169 

47.313 

0.170 

P  Aur 

54,403.8922 

0.878 

60.619 

0.289 

-98.087 

0.286 

P  Aur 

54,403.8982 

0.880 

61.415 

0.203 

-98.274 

0.196 

P  Aur 

54,404.8387 

0.117 

62.259 

0.292 

-99.681 

0.290 

P  Aur 

54,404.8542 

0.121 

60.690 

0.241 

-97.396 

0.236 

P  Aur 

54,404.8680 

0.125 

58.476 

0.253 

-95.754 

0.244 

P  Aur 

54,487.7015 

0.042 

86.798 

0.292 

-124.300 

0.274 

P  Aur 

54,487.7106 

0.044 

86.270 

0.237 

-124.136 

0.236 

P  Aur 

54,487.7191 

0.046 

86.305 

0.381 

-123.834 

0.374 

P  Aur 

54,487.7530 

0.055 

84.279 

0.550 

-122.404 

0.556 

P  Aur 

54,488.7643 

0.310 

-56.848 

0.217 

23.889 

0.209 

3 


The  Astronomical  Journal,  142:6  (llpp),  2011  July 


Table  2 

(Continued) 


Behr  et  al. 


Star  HJD  -  2,400,000 

Phase 

V\  (km  s  1 

)  Vi  Error  (km  s  1 ) 

V2  (km  s  ') 

Vi  Error  (km  s  1 ) 

fi  Aur 

54,488.7731 

0.312 

-58.307 

0.206 

24.076 

0.203 

fi  Aur 

54,488.7818 

0.315 

-59.895 

0.200 

26.191 

0.196 

fi  Aur 

54,491.6686 

0.044 

85.997 

0.288 

-124.270 

0.283 

fi  Aur 

54,491.6835 

0.047 

85.742 

0.159 

-123.753 

0.159 

fi  Aur 

54,491.7588 

0.066 

81.639 

0.350 

-119.413 

0.341 

fi  Aur 

54,491.7881 

0.074 

79.167 

0.413 

-116.398 

0.395 

Mizar 

54,488.9333 

0.660 

21.818 

0.361 

-37.796 

0.408 

Mizar 

54,488.9434 

0.661 

22.811 

0.980 

-37.749 

1.093 

Mizar 

54,488.9518 

0.661 

22.457 

0.312 

-37.301 

0.394 

Mizar 

54,488.9604 

0.662 

23.821 

0.556 

-37.302 

0.708 

Mizar 

54,492.0078 

0.810 

44.126 

0.247 

-57.208 

0.329 

Mizar 

54,492.0468 

0.812 

43.908 

0.354 

-57.704 

0.456 

Mizar 

54,548.8658 

0.578 

11.679 

0.482 

-26.290 

0.640 

Mizar 

54,549.8253 

0.625 

17.287 

0.174 

-32.974 

0.222 

Mizar 

54,576.7827 

0.937 

39.127 

0.417 

-52.615 

0.543 

Mizar 

54,578.7569 

0.034 

-74.826 

0.211 

59.789 

0.266 

Mizar 

54,579.7472 

0.082 

-81.159 

0.264 

66.865 

0.333 

Mizar 

54,579.7558 

0.082 

-81.404 

0.289 

66.645 

0.366 

Mizar 

54,602.7619 

0.202 

-52.427 

0.283 

37.381 

0.368 

Mizar 

54,603.6351 

0.245 

-42.830 

0.240 

27.669 

0.301 

Mizar 

54,604.6286 

0.293 

-33.460 

0.289 

19.960 

0.377 

Mizar 

54,634.6759 

0.756 

36.381 

0.306 

-49.740 

0.402 

Mizar 

54,636.6348 

0.851 

48.098 

0.184 

-61.666 

0.246 

Mizar 

54,637.6408 

0.900 

48.821 

0.213 

-62.867 

0.287 

k  Peg 

54,400.6227 

0.069 

26.262 

0.031 

2.603 

0.361 

k  Peg 

54,401.6788 

0.245 

-10.951 

0.035 

-0.613 

0.345 

k  Peg 

54,402.6158 

0.402 

-46.681 

0.055 

4.407 

0.415 

k  Peg 

54,402.6305 

0.405 

-47.035 

0.064 

3.843 

0.465 

k  Peg 

54,403.6430 

0.574 

-49.989 

0.039 

0.649 

0.309 

k  Peg 

54,404.5841 

0.732 

-16.839 

0.042 

-0.683 

0.416 

k  Peg 

54,404.5988 

0.734 

-16.173 

0.054 

0.520 

0.506 

k  Peg 

54,634.8531 

0.294 

-24.846 

0.029 

0.465 

0.317 

k  Peg 

54,635.8924 

0.468 

-55.088 

0.032 

-0.740 

0.284 

k  Peg 

54,636.8741 

0.632 

-41.896 

0.035 

6.552 

0.275 

k  Peg 

54,637.8349 

0.793 

-2.007 

0.034 

0.883 

0.333 

r]  Vir 

54,491.9611 

0.229 

11.665 

0.022 

-11.833 

0.067 

rj  Vir 

54,548.8113 

0.021 

-27.988 

0.034 

40.160 

0.101 

T]  Vir 

54,549.7519 

0.034 

-26.057 

0.030 

37.392 

0.090 

r)  Vir 

54,576.701 1 

0.409 

21.091 

0.035 

-25.264 

0.099 

rj  Vir 

54,577.7526 

0.424 

21.138 

0.035 

-25.293 

0.105 

rj  Vir 

54,579.7025 

0.451 

20.923 

0.039 

-25.080 

0.119 

r )  Vir 

54,603.7043 

0.785 

-6.225 

0.028 

10.650 

0.085 

r )  Vir 

54,604.6471 

0.799 

-8.312 

0.029 

13.633 

0.087 

r]  Vir 

54,605.6656 

0.813 

-10.717 

0.030 

16.643 

0.091 

r]  Vir 

54,634.6939 

0.217 

9.759 

0.033 

-10.524 

0.101 

ij  Vir 

54,636.6521 

0.244 

12.931 

0.026 

-14.812 

0.076 

Table  3 

Orbital  Parameters  and  Stellar  Mass  Estimates  for  i  Pegasi 

Parameter 

Fekel  &  Tomkin  (1983) 

This  Work  (Formal  RV  Errors) 

This  Work  (Scaled  RV  Errors) 

P  (days) 

10.213033  ±  0.000013 

adopted  from  F&T 

adopted  from  F&T 

T  (reduced  HJD) 

45320.1423 

54399.5296  ±  0.0003 

54399.5288  ±  0.0007 

e 

0.0  assumed 

0.0  assumed 

0.0  assumed 

K\  (km  s~*) 

48.1  ±  0.2 

48.380  ±  0.006 

48.380  ±  0.018 

K2  (km  s~ 1 ) 

77.9  ±  0.3 

77.637  ±  0.027 

77.638  ±  0.050 

V0  (km  s~!) 

-5.5  ±  0.2 

-4.245  ±  0.007 

-4.229  ±  0.015 

^obs 

32 

16 

16 

X2  primary 

142.20 

16.45 

X2  secondary 

64.52 

16.11 

(Jrv  primary  (km 

s-1) 

0.90 

0.056 

0.059 

(Jrv  secondary  (km  s  1 ) 

1.16 

0.158 

0.156 

Mi  (Mq  ) 

1.326  ±0.016a 

1.3239  ±  0.0018 

1.3239  ±  0.0025 

m2  (M0) 

0.819  ±  0.009a 

0.8250  ±  0.0010 

0.8250  ±  0.0013 

Note.  a  From  Boden  et  al.  (1999),  using  F&T  velocities  in  conjunction  with  spatial  interferometer  observations. 
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Vo  might  merely  be  a  result  of  different  RV  zero  points  between 
Fekel  &  Tomkin’s  observations  and  ours — unfortunately,  we  did 
not  make  any  RV  observations  of  i  Piscium,  the  RV  standard 
star  that  they  used  as  their  reference  spectrum. 

Although  the  rms  scatter  of  our  RV  points  around  the  best- 
fit  orbital  curves  is  small  (56  m  s"  1  and  158  m  s-1  for  the 
A  and  B  components,  respectively),  the  scatter  is  larger  than 
would  be  expected  from  the  error  bars  on  each  individual  RV 
measurement,  and  is  significantly  above  the  instrumental  RV 
error  floor  of  ~10  m  s^1  that  we  determined  in  Behr  et  al. 
(2009).  To  account  for  this  discrepancy,  we  multiply  the  RV 
error  bars  by  3.11  for  the  primary  and  1.83  for  the  secondary, 
such  that  the  mean  per-measurement  error  bar  matches  the  rms 
deviation  ctrv  for  each  stellar  component.  The  orbital  fits  are 
then  recalculated  using  these  scaled  error  bars,  and  the  resulting 
orbital  parameters  are  listed  in  the  rightmost  column  of  Table  3. 
This  same  procedure  is  applied  to  all  subsequent  binary  systems 
as  well. 

The  additional  RV  variability,  if  real,  may  be  a  result  of 
stellar  activity  on  both  the  primary  and  secondary,  driven 
by  tidal  interactions  between  the  two  stars.  Fekel  &  Tomkin 
measure  v  sin  /  =  7  ±  2  km  s"1  for  the  primary,  very  close 
to  an  estimated  synchronous  rate  of  6.5  km  s_1,  and  v  sin  i  = 
9±3  km  s“ 1  for  the  secondary,  which  is  well  above  the  estimated 
synchronous  rate  of  4.5  km  s"1 .  Gray  (1984)  finds  v  sin  i  values 
of  6.5  ±  0.3  km  s"1  (primary)  and  5  ±  1  km  s"1  (secondary), 
suggesting  that  the  system  is  synchronized.  Our  preliminary 
analysis  of  line  broadening  indicates  v  sin/  =  7.6  and  7.2  for 
the  primary  and  secondary,  respectively.  If  the  secondary  is 
indeed  spinning  more  rapidly  than  the  synchronous  rate,  then 
above-average  surface  activity  could  result,  which  would  add 
significant  astrophysical  RV  “jitter”  to  our  measurements.  A 
synchronously  rotating  component  would  be  less  susceptible 
to  tidal  effects,  but  activity  might  still  be  enhanced  by  the 
proximity  of  a  massive  companion.  However,  Konacki  et  al. 
(2009)  measured  RVs  of  i  Peg  with  three  different  spectrographs 
and  found  no  jitter  greater  than  17ms"1  (primary)  and  85  ms"1 
(secondary),  suggesting  that  the  jitter  observed  by  dFTS2  was 
instrumental  rather  than  astrophysical. 

Boden  et  al.  (1999)  measured  an  inclination  angle  for  this 
system  of  i  —  95?67  ±  0122  (based  on  their  primary  data  set). 
Using  this  value  along  with  our  orbital  parameters,  with  the 
fundamental  parameters  recommended  by  Torres  et  al.  (2010), 
we  derive  stellar  masses  of  M\  =  1.3241  ±  0.0018  Mq  and 
Mi  —  0.825 1  ±0.0010  Mq,  which  represent  relative  (statistical) 
errors  of  0. 14%  and  0. 12%,  respectively.  With  scaled  error  bars, 
the  mass  estimates  are  the  same,  albeit  with  larger  error  bars, 
for  relative  uncertainties  of  0.19%  and  0.15%.  These  values 
agree  reasonably  well  with  the  calculations  of  Boden  et  al. 
who  used  Fekel  &  Tomkin’s  K\  and  Ki  values  to  determine 
M,  =  1.326  ±  0.016  Mq  and  M2  =  0.819  ±0.009  M0.  For  our 
mass  estimates,  the  largest  component  of  the  error  budget  is  due 
to  the  uncertainty  in  /,  although  the  uncertainties  in  the  K  values 
are  also  significant  contributors. 

3.2.  co  Draconis 

The  spectroscopic  orbit  of  the  co  Draconis  system  (HR  6596, 
HD  160922,  and  HIP  86201)  was  measured  by  Mayor  &  Mazeh 
(1987)  and  more  recently  by  Fekel  et  al.  (2009).  Their  derived 
orbital  parameters  are  shown  in  Table  4,  along  with  our  values. 
Our  K  velocities  agree  closely  with  those  of  Fekel  et  al.  We  find 
a  small  but  non-zero  eccentricity  for  the  orbits,  0.0023  ±  0.0002. 
Fekel  et  al.  derived  e  —  0.0027  ±  0.0008  for  the  primary,  with 


Figure  2.  RV  measurements  of  the  SB2  system  co  Draconis. 


coi  =  40: 1  ±  1 7:8  fF.  C.  Fekel  2009,  private  communication  ),  but 
their  e  and  co  values  for  the  secondary  did  not  agree  with  those 
of  the  primary,  so  they  adopted  a  circular  orbit  for  the  system. 
(Our  co\  =  137.86  ±  13.48,  as  described  below.)  Despite  the 
similarity  between  our  e  value  and  their  non-zero  e  value,  the 
measurements  of  co\  are  substantially  different,  so  we  cannot 
plausibly  claim  that  an  orbital  eccentricity  has  been  clearly 
detected. 

As  an  additional  test  of  the  non-zero  eccentricity,  we  follow 
the  Fekel  et  al.  procedure  of  fitting  orbits  to  the  A  and  B 
components  separately  and  comparing  the  derived  values  for 
co\  and  ®2,  which  should  differ  by  180°.  Using  the  formal 
error  bar  data,  we  calculated  e\  —  0.0020  ±  0.0003  with 
co\  =  74°  16  ±  12?74,  and  e 2  =  0.0027  ±  0.0005  with  C02  = 
243?  1 4  ±  1 3?05.  The  eccentricity  values  agree  reasonably  well, 
and  the  co  angles  differ  by  ~169°,  which  is  within  1  a  of  180°. 
However,  this  co\  value  does  not  agree  with  the  o>\  value  from 
the  combined  fit,  which  is  puzzling.  This  discrepancy  may  be 
related  to  the  apparent  systematic  trends  in  the  secondary  RV 
residuals  which  are  evident  in  Figure  2.  Further  observations 
with  better  phase  coverage  will  be  required  to  validate  or  refute 
our  measured  eccentricity  for  co  Dra. 

As  with  i  Peg,  our  value  for  Vo  differs  from  the  prior  work  by 
a  statistically  significant  amount,  although  the  magnitude  of  the 
difference  is  not  as  large.  Differences  in  the  RV  zero  point  are 
the  most  likely  explanation. 

No  visual  orbit  or  estimate  of  the  inclination  angle  i  has 
been  determined  for  this  binary,  despite  efforts  to  resolve  it 
using  speckle  interferometry  (Isobe  1991;  Miura  et  al.  1995). 
We  are  therefore  unable  to  calculate  the  true  masses  of  the 
stellar  components.  We  find  M \  sin3  i  —  0.16054±0.0001 1  Me 
and  Mi  sin3  i  =  0.13030  ±  0.00007  Mq,  in  moderately  good 
agreement  with  Fekel  et  al.  (With  the  scaled  RV  error  bars, 
our  mass  estimates  are  virtually  unchanged,  with  error  bars 
approximately  three  times  larger.)  We  hope  that  long-baseline 
interferometers  will  soon  be  able  to  resolve  the  astrometric  orbit 
of  this  system  and  determine  the  inclination  angle. 


5 


The  Astronomical  Journal,  142:6  (llpp),  2011  July 


Behr  et  al. 


Table  4 

Orbital  Parameters  and  Stellar  Mass  Estimates  for  u>  Draconis 


Parameter 

Mayor  &  Mazeh  (1987) 

Fekel  et  al.  (2009) 

This  Work  (Formal  RV  Errors) 

This  Work  (Scaled  RV  Errors) 

P  (days) 

5.279799  ±  0.000003 

5.2798088  ±  0.0000083 

adopted  from  Fekel  et  al. 

adopted  from  Fekel  et  al. 

T  (reduced  HJD) 

44698.273  ±  0.005 

53980.1606  ±  0.0006 

54547.1347  ±  0.0753 

54547.1180  ±  0.1974 

e 

0 

0.0  assumed 

0.0023  ±  0.0002“ 

0.0023  ±  0.0006 

tol  (deg) 

139.01  ±  5.14 

137.86  ±  13.48 

K\  (km  s~') 

35.8  ±  0.3 

36.326  ±  0.029 

36.293  ±  0.008 

36.292  ±  0.020 

Ko  (km  s~ 1 ) 

45.2  ±  0.3 

44.699  ±  0.039 

44.717  ±  0.014 

44.718  ±  0.038 

Vo  (km  s~‘) 

-14.1  ±  0.2 

-13.975  ±  0.018 

-13.497  ±  0.006 

-13.501  ±  0.016 

Mobs 

27 

82 

16 

16 

X2  primary 

92.10 

15.42 

X2  secondary 

138.65 

19.42 

<trv  primary  (km  s-1) 

0.19  (unit  weight) 

0.061 

0.059 

ctrv  secondary  (km  s-1) 

0.126 

0.126 

Mi  sin3  i  (Mq) 

0.163  ±  0.003 

0.16090  ±  0.00032 

0.16054  ±  0.00011 

0.16054  ±  0.00030 

M2  sin3  i  (Mq) 

0.129  ±  0.002 

0.13076  ±  0.00024 

0.13030  ±  0.00007 

0.13029  ±  0.00018 

Note.  a  But  see  the  text  regarding  the  validity  of  this  non-zero  eccentricity. 

Table  5 

Orbital  Parameters  and  Stellar  Mass  Estimates  for  12  Bootis 

Parameter 

Boden  et  al.  (2005)a 

Tomkin  &  Fekel  (2006) 

This  Work  (Formal  RV  Errors) 

This  Work  (Scaled  RV  Errors) 

P  (days) 

9.6045492  ±  0.0000076 

9.6045529  ±  0.0000048 

adopted  from  T&F 

adopted  from  T&F 

T  (reduced  HJD) 

51237.7729  ±  0.0051 

52400.4292  ±  0.0035 

54542.2431  ±  0.0031 

54542.2424  ±  0.0042 

e 

0.19233  ±  0.00086 

0.19268  ±  0.00042 

0.1928  ±  0.0003 

0.1928  ±  0.0004 

a>l  (deg) 

286.67  ±  0.19 

286.87  ±  0.14 

286.79  ±  0.12 

286.78  ±  0.17 

A'i  (km  s~*) 

67.302  ±  0.087 

67.286  ±  0.037 

67.107  ±  0.035 

67.113  ±  0.047 

Ko  (km  s~ 1 ) 

69.36  ±  0.10 

69.30  ±  0.05 

69.110  ±  0.037 

69.102  ±  0.054 

V0  (km  s~') 

9.551  ±  0.051 

9.578  ±  0.022 

10.040  ±  0.018 

10.046  ±  0.025 

M>bs 

49 

24 

13 

13 

X2  primary 

-49.0 

17.87 

10.16 

X2  secondary 

-49.0 

30.79 

14.02 

(Try  primary  (km  s-1) 

0.47 

0.11 

0.082 

0.077 

<7ry  secondary  (km  s-1) 

0.54 

0.115 

0.119 

Ml  (Mq) 

1.4160  ±  0.0049 

1.416  ±  0.003 

1.4013  ±  0.0025 

1.4011  ±  0.0031 

Mo  (Mq) 

1.3740  ±  0.0045 

1.375  ±  0.002 

1.3607  ±  0.0024 

1.3608  ±  0.0028 

Note.  a  Combined  fit  to  RV  and  astrometric  data. 


3.3.  12  Bodtis 

The  spectroscopic  binary  12  Bootis  (HR  5304,  HD  123999, 
and  HIP  69226)  has  received  recent  attention  from  both  Boden 
et  al.  (2005),  who  combined  spectroscopic  and  astrometric  data, 
and  Tomkin  &  Fekel  (2006),  who  performed  a  high-precision 
spectroscopy-only  assessment  of  the  orbit.  Orbital  parameters 
are  shown  in  Table  5,  and  our  RV  data  are  plotted  in  Figure  3. 
The  derived  quantities  for  e  and  a>\  are  in  excellent  agreement 
among  all  three  studies.  Our  K\  and  Ko  values,  on  the  other 
hand,  are  smaller  than  those  of  Boden  et  al.  and  Tomkin  & 
Fekel  by  several  standard  deviations,  and  our  derived  systemic 
velocity  is  different  as  well.  The  discrepancy  in  Vo  may  simply 
be  ascribed  to  a  different  RV  zero  point,  but  the  difference  in 
K\/2  deserves  further  scrutiny.  Due  to  the  premature  conclusion 
of  our  observing  program,  our  RV  data  do  not  fully  cover  the 
region  of  maximum  absolute  velocities  around  phase  =  0.15, 
so  the  velocity  amplitudes  are  not  as  reliably  constrained  as 
they  might  be.  When  dFTS  observations  resume,  12  Boo  will 
be  one  of  our  highest  priority  targets,  so  that  this  issue  can  be 
addressed. 

Given  smaller  K  amplitudes  than  prior  publications,  we  derive 
smaller  masses  as  well.  Using  i  =  107°990  ±  0°077  from 
Boden  et  al.,  we  determine  M\  =  1.4013  ±  0.0025  Mq  and 
Mo  =  1.3607  ±  0.0024  Mq,  for  a  relative  statistical  uncertainty 
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phase 

Figure  3.  RV  measurements  of  the  SB2  system  12  Bootis. 
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Figure  4.  RV  measurements  of  the  SB2  system  VI 143  Cygni. 


Figure  5.  RV  measurements  of  the  SB2  system  (3  Aurigae. 


of  0.18%.  These  values  are  several  a  smaller  than  the  masses 
derived  by  Boden  et  al.  and  Tomkin  &  Fekel. 


3.4.  VI 143  Cygni 

The  eclipsing  SB2  system  of  VI 143  Cygni  (HR  7484, 
HD  185912,  and  HIP  96620)  was  previously  analyzed  by 
Andersen  et  al.  (1987).  Their  orbital  parameters  are  compared 
to  ours  in  Table  6,  and  our  RV  data  are  plotted  in  Figure  4. 
There  is  a  broad  agreement  between  the  two  sets  of  orbital 
elements,  although  the  error  bars  that  we  derived  for  K\  and 
/G  are  relatively  large,  partly  because  of  the  small  number  of 
observations,  partly  because  the  individual  RV  measurements 
had  larger  error  bars  due  to  larger  rotational  broadening  of  the 
absorption  lines  and  lower  signal-to-noise  ratio.  Of  particular 
interest  is  the  comparison  of  the  periastron  angle  co\.  From 
precise  photometric  timing  of  the  system’s  eclipses,  Gimenez 
&  Margrave  (1985)  detected  apsidal  motion  (precession  of  the 
periastron  point)  with  a  period  of  10,750  years,  which  would 
imply  a  change  in  the  value  of  u>\  of  0.76  during  the  ~22.6  years 
that  elapsed  between  their  last  observations  (1985  October)  and 
our  first  observations  (2008  May).  The  actual  measured  change 
in  a>\  is  +0°39  ±  (K25.  The  predicted  periastron  precession 
is  therefore  not  ruled  out,  but  is  not  solidly  confirmed  either. 
With  more  extensive  observations  of  VI 143  Cyg,  we  hope  to 
place  more  useful  constraints  on  the  magnitudes  of  the  classical 
gravitational  quadrupole  and  general  relativity  effects  which 
cause  the  precession. 

Andersen  et  al.  (1987)  adopt  an  inclination  angle  of  87.0  ±  1° 
based  upon  eclipse  photometry  by  Wood  (1971),  Popper  & 
Etzel  (1981),  and  van  Hamme  &  Wilson  (1984).  Using  that 
same  value  for  i,  we  estimate  M\  =  1.3815  ±  0.01 14  MG 
and  M2  =  1.3451  ±  0.0100  MG.  The  error  bars  on  K\  and 
K2  dominate  the  error  budget  for  the  masses,  so  further  high- 
accuracy  spectroscopic  observations  of  this  system  are  clearly 
called  for. 


3.5.  f3  Aurigae 

ft  Aurigae  (HR  2088,  HD  40183,  and  HIP  28360)  is  another 
eclipsing  double-lined  binary,  consisting  of  two  A2  subgiants  in 
a  close  four-day  orbit.  Smith  (1948)  measured  the  RV  curves  of 
both  components,  and  Pourbaix  (2000)  reanalyzed  these  data  in 
conjunction  with  interferometric  astrometry  data  from  Hummel 
et  al.  (1995)  to  refine  the  orbital  parameters.  Our  RV  curve  is 
displayed  in  Figure  5.  The  phase  coverage  was  insufficient  to 
constrain  e  or  coy  so  we  assumed  a  circular  orbit.  Our  derived 
parameters  differ  from  the  prior  two  analyses  (Table  7),  with  K\ 
and  K2  semi-amplitude  values  intermediate  between  those  of 
Smith  and  those  of  Pourbaix.  As  with  V 1 143  Cyg,  this  system 
exhibits  a  rotational  line  broadening  of  30-40  km  s-1,  which 
increases  the  uncertainty  of  each  RV  measurement  and  thus  the 
derived  orbital  parameters. 

Adopting  i  —  76:0  ±  0.4  from  Hummel  et  al.  (1995),  we 
calculate  Mi  =  2.3885  ±  0.0134  MG  and  M2  =  2.3270  ± 
0.0130  Mg.  (Pourbaix  uses  i  =  75°0  ±  0°73;  the  source  of 
this  value  is  unclear.)  The  majority  of  the  0.54%  relative  error 
in  our  mass  values  is  due  to  the  inclination  angle  uncertainty, 
so  this  system  would  be  a  prime  follow-up  target  for  further 
long-baseline  spatial  interferometry. 

3.6.  Mizar  A 

The  brighter  component  of  the  visual  binary  Mizar  (HR 
5054,  HD  116656,  and  HIP  65378)  is  itself  a  spectroscopic 
binary,  with  two  early-A  dwarfs  in  an  elliptical  20  day  orbit. 
Table  8  lists  the  orbital  parameters  measured  by  Fehrenbach 
&  Prevot  (1961)  and  the  subsequent  revisions  computed  by 
Pourbaix  (2000)  with  astrometry  data  from  Hummel  et  al. 
(1998),  along  with  the  values  computed  from  our  RV  data 
(Figure  6).  We  find  smaller  velocity  amplitudes  and  a  slightly 
larger  eccentricity  for  this  binary  system  than  prior  researchers. 
Like  the  prior  two  targets,  Mizar  A’s  component  spectra  are 
moderately  rotationally  broadened,  reducing  the  quality  of  the 
RV  measurements. 
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Table  6 

Orbital  Parameters  and  Stellar  Mass  Estimates  for  VI 143  Cygni 


Parameter 

Andersen  et  al.  (1987)  This  Work  (Formal  RV  Errors) 

This  Work  (Scaled  RV  Errors) 

P  (days) 

7.64075217  ± 

0.0000005 1  adopted  from  Andersen 

adopted  from  Andersen 

T  (reduced  HJD) 

42212.76652  ± 

0.00015  54598.1853  ±  0.0027 

54598.1835  ±  0.0088 

e 

0.540  ± 

0.003 

0.5469  ±  0.0010 

0.5484  ±  0.0032 

nil  (deg) 

48.6  ± 

0.02 

48.99  ±  0.25 

48.84  ±  0.83 

K\  (km  s-1) 

88.20  ± 

0.20 

88.867  ±  0.248 

89.055  ±  0.836 

K2  (km  s' 1 ) 

91.10  ± 

0.40 

91.267  ±  0.311 

91.508  ±  0.941 

V0  (km  s'1) 

-16.5  ± 

0.7 

-16.505  ±  0.074 

-16.461  ±  0.235 

A^obs 

62 

7 

7 

X2  primary 

63.01 

5.98 

X2  secondary 

46.49 

5.04 

(Jrv  primary  (km  s_1) 

1.1 

0.734 

0.789 

(Try  secondary  (km  s_ 

) 

2.2 

1.012 

0.951 

Mi  (M0) 

1.391  ± 

0.016 

1.3815  ±  0.0114 

1.3868  ±  0.0334 

M2  (Mo) 

1.347  ± 

0.013 

1.3451  ±  0.0100 

1.3496  ±  0.0308 

Table  7 

Orbital  Parameters  and  Stellar  Mass  Estimates  for  f$  Aurigae 

Parameter 

Smith  (1948) 

Pourbaix  (2000)a 

This  Work  (Formal  RV  Errors) 

This  Work  (Scaled  RV  Errors) 

P  (days) 

3.9600421  ±  0.0000013 

3.96004  ±  0.00000267 

adopted  from  Pourbaix 

adopted  from  Pourbaix 

T  (reduced  HJD) 

31076.719 

43915.7 

54539.0162  ±  0.0003 

54537.0362  ±  0.0004 

e 

0.0 

2.75266  x  10'6  ±  0.007 

0.0  assumed 

0.0  assumed 

O),  (deg) 

0.0 

139.043  ±  360.0 

K\  (km  s'1) 

107.46  ±  0.39 

110.246  ±  1 

108.053  ±  0.072 

108.053  ±  0.099 

K2  (km  s' 1 ) 

111.49  ±  0.37 

110.52  ±  2.1 

110.911  ±  0.071 

110.911  ±  0.098 

Vo  (km  s'1) 

-17.06  ±  0.27 

-15.7536  ±  0.62 

-17.552  ±  0.037 

-17.552  ±  0.052 

Mobs 

21 

21 

20 

20 

X2  primary 

54.19 

28.30 

X2  secondary 

40.83 

21.15 

ctrv  primary  (km  s-1) 

2.740 

0.363 

0.363 

(Try  secondary  (km  s“ 

) 

6.369 

0.358 

0.358 

Mi  (M0) 

2.4  ±  0.1 

2.3885  ±  0.0129 

2.3885  ±  0.0134 

M2  (Mo) 

2.44  ±  0.073 

2.3270  ±  0.0126 

2.3270  ±  0.0130 

Note.  a  These  orbital  parameters  are  listed  in  the  downloadable  data  table  at  the  SB9  Web  site  (http://sb9.astro.ulb.ac 

be/mainform.cgi)  but  cannot  be 

accessed  directly  from  the  Web  interface. 

Table  8 

Orbital  Parameters  and  Stellar  Mass  Estimates  for  Mizar  A 

Parameter 

Fehrenbach  &  Prevot  (1961)  Pourbaix  (2000)a 

This  Work  (Formal  RV  Errors) 

This  Work  (Scaled  RV  Errors) 

P  (days) 

20.5386 

20.5385  ±  0.00013514 

adopted  from  Pourbaix 

adopted  from  Pourbaix 

T  (reduced  HJD) 

36997.212  ±  0.022 

38085.7  ±  0.0269224 

54536.9882  ±  0.0068 

54536.9904  ±  0.0106 

e 

0.537  ±  0.004 

0.529404  ±  0.0052 

0.5415  ±  0.0010 

0.5415  ±  0.0016 

0)1  (deg) 

104.16  ±  1.15 

105.5  ±  0.79 

105.21  ±  0.14 

105.27  ±  0.23 

Ki  (km  s'1) 

68.80  ±  0.79 

67.2586  ±  0.96 

66.479  ±  0.095 

66.478  ±  0.153 

K2  (km  s' 1 ) 

67.60  ±  0.91 

69.1796  ±  0.77 

66.012  ±  0.118 

66.019  ±  0.177 

Vo  (km  s'1) 

-5.64  ±  0.15 

-6.3077  ±  0.38 

-7.342  ±  0.052 

-7.309  ±  0.081 

Mobs 

15 

15 

18 

18 

X2  primary 

82.55 

33.70 

X2  secondary 

81.44 

34.01 

(Try  primary  (km  s-1) 

1.87 

1.88094 

0.556 

0.566 

(Try  secondary  (km  s“ 

) 

1.32 

2.39922 

0.641 

0.638 

Mi  (Mo) 

2.5  ±  0.11 

2.2224  ±  0.0221 

2.2228  ±  0.0250 

M2  (Mo) 

2.5  ±  0.12 

2.2381  ±  0.0219 

2.2383  ±  0.0246 

Note.  a  With  some  values  from  the  SB9  catalog  (Pourbaix  et  al.  2004). 


According  to  Hummel  et  al.,  i  =  60.5  ±  053  for  this  binary. 
Using  this  value,  we  determine  that  Mi  =  2.2224  ±  0.0221  MG 
and  M2  =  2.2381  ±  0.0219  M0,  for  a  relative  mass  error 
of  1.00%  and  0.98%,  respectively.  The  uncertainty  in  i  is 
responsible  for  most  of  the  mass  error;  if  a,  could  be  reduced  to 
0505  and  the  errors  on  K\  and  Ko  could  be  cut  in  half,  then  the 
mass  uncertainty  would  drop  below  0.25%. 


4.  RESULTS  FOR  TRIPLE  SYSTEMS 
4. 1.  k  Pegasi 

The  k  Pegasi  system  (HR  83 15,  HD  206901,  and  HIP  107354) 
is  a  hierarchical  triple,  with  two  bright  components  (A  and  B) 
in  an  11.5  year  orbit,  and  a  fainter  unseen  component  in  a 
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Table  9 

Orbital  Parameters  for  the  Single-lined  B  Component  of  k  Pegasi 


Parameter 

Mayor  &  Mazeh  (1987) 

Hajian  et  al.  (2007) 

Muterspaugh  et  al.  (2008) 

This  Work  (2007  Oct) 

This  Work  (2008  June) 

P  (days) 

5.97164  ±  0.00006 

adopted  from  M&M 

5.9714971  ±  0.0000013 

adopted  from  M&M 

adopted  from  M&M 

T  (reduced  HJD) 

44801.589  ±  0.015 

53681.86  ±  0.04 

52402.22  ±  0.10 

54400.2125  ±  0.0005 

54633.0989  ±  0.0005 

e 

0 

adopted  from  M&M 

0.0073  ±  0.0013 

adopted  from  M&M 

adopted  from  M&M 

K\  (km  s_1) 

42.1  ±  0.3 

41.572  ±  0.257 

42.301  ±  0.025 

42.527  ±  0.037 

V0  (km  s_  1 ) 

-0.8  ±  0.2 

-9.40  ±  0.22 

-12.135  ±  0.016 

-13.333  ±  0.022 

VA  (km  s_1) 

1.352  ±  0.851 

2.019  ±  1.876 

Nobs 

30 

9 

30 

7 

4 

c7rv  B  (km  s_1) 

1.1 

0.990 

0.035 

0.063 

0.093 

<trv  A  (km  s"1) 

0.250 

2.084 

3.249 

secondory 

5  i 

5 

^  i 

*  M  5  5 

. 5  5  I 

- 

0.0 
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Figure  6.  RV  measurements  of  the  SB2  system  Mizar. 


six-day  orbit  around  the  B  component.  The  canonical  published 
orbit  for  this  system  comes  from  Mayor  &  Mazeh  (1987),  with 
more  recent  observations  by  Konacki  (2005)  (with  orbit  analyses 
published  in  Muterspaugh  et  al.  2006  and  Muterspaugh  et  al. 
2008)  and  our  dFTSl  prototype  (Hajian  et  al.  2007).  Table  9 
displays  the  orbital  parameters  as  measured  by  Mayor  &  Mazeh, 
dFTSl,  and  dFTS2.  Following  Mayor  &  Mazeh,  we  assume 
that  e  —  0  for  the  short-period  orbit.  The  dFTS2  observations 
of  this  system  were  made  during  two  observing  runs  separated 
by  ~0.64  years,  so  we  expect  the  Vo  value  of  the  short-period 
binary  to  change  due  to  the  long-period  orbit.  To  account  for 
this  change,  we  treated  the  data  from  the  two  observing  runs 
completely  separately,  and  the  results  are  given  in  two  separate 
columns  in  the  table.  For  the  RV  plot  in  Figure  7,  we  shifted 
Bb  component’s  RV  values  from  the  first  observing  run  by 
—  1.198  km  s-1  so  that  the  two  different  epochs  would  share 
the  same  Vo. 

With  only  four  RV  measurements  in  the  second  observing  run, 
the  value  for  the  K\  amplitude  should  be  considered  provisional, 
but  the  change  in  V()  is  quite  clear.  We  see  an  even  larger 
change  in  Vo  as  compared  to  the  mid-1981  observations  of 
Mayor  &  Mazeh,  although  differences  in  RV  zero  point  must 
be  considered.  Unfortunately,  our  prior  dFTS  1  observations  did 


;  A  component  n 


§  0'~ . u- 

■o  ;  q  D  $  s 

S  “2  “  n  n  5 
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Figure  7.  RV  measurements  of  the  double-lined  triple  system  k  Pegasi.  Filled 
symbols  denote  the  primary  of  the  short-period  pair  (Ba),  while  open  symbols 
denote  the  “A”  component  of  the  long-period  orbit. 

not  yield  a  value  of  Vo,  because  the  template  spectra  were  not 
referenced  to  an  absolute  wavelength  standard. 

RV  measurements  of  the  A  component  are  not  as  precise  as 
those  of  the  B  component,  because  its  lines  are  much  broader: 
we  estimate  usini  =  47.9  km  s_1  for  A  and  7.1  km  s'1  for 
B.  Even  taking  this  fact  into  account,  however,  we  find  a  much 
larger  scatter  of  our  RV  measurements  for  A  than  expected  from 
X2  statistics.  One  possible  explanation  for  this  discrepancy  is 
that  the  A  component  is  also  a  close  binary,  as  proposed  by 
Beardsley  &  King  (1976).  We  phased  our  RV  data  to  their 
claimed  4.77  day  period,  but  did  not  find  any  coherent  cyclic 
pattern  in  the  A  component  velocities.  (Muterspaugh  et  al.  2006 
see  no  evidence  for  a  fourth  component  either.)  The  binarity  of 
A  is  a  possibility  that  we  might  explore  with  future  data,  but  for 
the  time  being,  this  hypothesis  is  not  supported. 

4.2.  i]  Virginis 

The  ii  Virginis  triple  system  (HR  4689,  HD  107259,  and 
HIP  60129)  was  studied  extensively  by  Hartkopf  et  al.  (1992) 
who  combined  spectroscopy  and  speckle  interferometry  to  de¬ 
termine  the  orbits  of  both  the  short-period  (72  day)  pair  and 
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Table  10 

Orbital  Parameters  and  Stellar  Mass  Estimates  for  t;  Virginis  Aa-Ab 


Parameter 

Hartkopf  et  al.  (1992) 

Hummel  et  al.  (2003) 

This  Work  (Formal  RV  Errors) 

This  Work  (Scaled  RV  Errors) 

P  (days) 

71.7919  ±  0.0009 

71.7916  ±  0.0006 

adopted  from  Hartkopf  et  al. 

adopted  from  Hartkopf  et  al. 

T  (reduced  HJD) 

47583.98  ±  0.25 

52321.4  ±  0.3 

54403.7295  ±  0.0938 

54403.6116  ±  0.3406 

e 

0.272  ±  0.009  (Aa) 
0.258  ±  0.012  (Ab) 

0.244  ±  0.007 

0.2518  ±  0.0011 

0.2519  ±  0.0040 

0)l  (deg) 

200.9  ±  1.5 

196.9  ±  1.8 

197.96  ±  0.48 

197.21  ±  1.74 

Ki  (km  s~') 

26.67  ±  0.20 

Same  as  Hartkopf  et  al. 

26.532  ±  0.054 

26.606  ±  0.198 

Sa  (km  s~ 1 ) 

35.58  ±  0.31 

Same  as  Hartkopf  et  al. 

35.128  ±  0.081 

35.236  ±  0.273 

Vo  (km  s_l) 

5.24  ±  0.19  (Aa) 
4.85  ±  0.32  (Ab) 

4.9  ±  0.2 

1.055  ±  0.009 

1.118  ±  0.033 

^obs 

50 

Same  as  Hartkopf  et  al. 

11 

11 

X2  primary 

292.62 

25.75 

X2  secondary 

263.24 

28.37 

(Try  primary  (km  s_1) 

1.96 

Same  as  Hartkopf  et  al. 

0.129 

0.150 

ctry  secondary  (km  s-1) 

4.12 

Same  as  Hartkopf  et  al. 

0.231 

0.218 

Mi  (M0) 

2.34  ±  0.2 

2.68  ±  0.15 

2.4818  ±  0.1158 

2.5039  ±  0.1246 

Mi  (Mo) 

1.95  ±  0.2a 

2.04  ±  0.10 

1.8745  ±  0.0874 

1.8907  ±  0.0932 

Note.  a  Assumed. 


Figure  8.  RV  measurements  of  the  double-lined  triple  system  ij  Virginis. 


the  long-period  (13.1  year)  grouping.  Hummel  et  al.  (2003) 
made  additional  observations  of  this  system  with  the  NPOI  in¬ 
terferometer,  refining  the  orbital  parameters  and  determining  the 
inclination  angle  of  the  close  binary  orbit.  Table  10  compares 
their  parameters  for  the  short-period  Aa-Ab  pair  to  our  deriva¬ 
tion.  The  parameters  are  in  general  agreement,  except  that  the 
change  in  Vo  is  larger  than  previously  seen,  and  is  likely  due  to 
the  gravitational  influence  of  the  third  component  of  the  system. 
Figure  8  shows  the  RV  curves  for  the  close  binary,  as  measured 
by  dFTS2.  Note  that  the  residuals  for  the  secondary  RV  points 
all  lie  above  the  dotted  line  denoting  zero  residual.  This  offset 
may  indicate  a  significant  difference  in  the  gravitational  redshift 
or  convective  blueshift  between  the  two  component  stars,  or  it 
may  be  an  effect  of  the  third  component.  Interestingly,  Hartkopf 
et  al.  see  a  similar  effect,  but  with  the  opposite  sign,  in  that  the 


Vo  that  they  derive  from  the  Aa  RV  curve  is  ~0.4  km  s-1  larger 
than  Vo  from  the  Ab  component. 

Hartkopf  et  al.  made  a  tentative  spectroscopic  detection  of 
a  blended  Mgn  4481  feature  from  the  faint  tertiary  compo¬ 
nent,  which  appeared  to  be  rotationally  broadened  by  about 
160  km  s-1.  Our  instrument  bandpass  does  not  include  this 
line,  so  we  were  unable  to  verify  their  detection,  but  we  per¬ 
formed  a  crude  three-dimensional  cross-correlation  using  two 
narrow-lined  template  spectra  plus  abroad-lined  A2V  template. 
We  found  no  evidence  of  a  consistent  RV  solution  for  the  third 
component. 

According  to  the  observations  of  Hummel  et  al.  (2003),  the 
inclination  angle  of  r)  Vir  Aa-Ab  is  45°5  ±  0.9.  Combining 
this  number  with  our  orbital  parameters,  we  determine  that 
Mi  =  2.4818  ±  0.1158  M0  and  M2  =  1.8745  ±  0.0874  M0. 
The  inclination  uncertainty  is  the  dominant  contributor  to  the 
error  budget,  although  the  velocimetry  results  can  certainly  be 
improved  with  more  data  points. 

5.  SUMMARY 

Our  results  demonstrate  that  dFTS  technology  is  well  suited  to 
high-accuracy  RV  measurements  of  double-lined  spectroscopic 
systems.  We  have  determined  the  orbital  parameters  of  six 
binary  systems,  matching  or  improving  the  published  values 
for  the  masses  of  the  component  stars.  We  also  observed  two 
double-lined  triple  systems,  providing  some  constraints  on  the 
nature  of  their  stars. 

For  our  future  observational  programs  for  spectroscopic  bi¬ 
nary  stars,  we  are  motivated  by  an  assortment  of  specific  scien¬ 
tific  goals  for  which  the  capabilities  of  a  dFTS  are  particularly 
applicable. 

1 .  The  most  immediate  goal  is  to  continue  to  improve  the  ac¬ 
curacy  of  orbital  parameters  of  binary  systems,  particularly 
the  K  amplitudes,  and  thus  measure  stellar  masses  more  ac¬ 
curately.  These  advancements  in  spectroscopic  capabilities 
must  proceed  in  parallel  with  better  astrometric  measure¬ 
ments,  as  determined  by  current  and  future  long-baseline 
spatial  interferometers. 

2.  With  high-accuracy  RV  measurements  spanning  longer 
periods  of  time,  we  will  be  able  to  detect  and  quantify 
secular  changes  in  binaries’  systemic  velocities  (Vo)  due 
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to  tertiary  companions.  As  discussed  by  Tokovinin  et  al. 
(2006),  the  presence  of  a  tertiary  companion  has  significant 
implications  for  the  formation  of  close  binaries. 

3.  Observations  of  near-circular  binary  orbits  will  confirm 
or  refute  small  non-zero  eccentricities,  thus  providing 
observational  validation  for  theories  of  tidal  circularization 
and  the  influence  of  external  gravitational  perturbations 
such  as  Kozai  resonances. 

4.  For  highly  elliptical  systems  like  VI 143  Cyg,  long-term 
observing  programs  can  measure  changes  in  periastron 
angle  to  test  theories  of  apsidal  precession  due  to  classical 
and  relativistic  effects. 

5.  Because  the  instrumental  profile  of  a  dFTS  is  easy  to 
calculate  a  priori  from  the  delay  sampling  function,  we 
can  measure  spectral  line  broadening  very  accurately,  e.g., 
to  determine  the  projected  rotational  velocities  of  stellar 
components  and  thus  shed  light  on  tidal  spin-up/spin-down 
mechanisms. 

This  research  was  funded  in  part  by  ARH’s  NSERC  Discovery 
Grant. 
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tion  of  dFTS2,  and  we  thank  telescope  operators  Geno  Bechetti, 
Dennis  Means,  and  Peter  Milne  for  their  expertise  in  operating 
the  telescope  on  our  behalf.  We  also  express  our  appreciation 
to  the  Director  of  the  Steward  Observatory  for  granting  us  tele¬ 
scope  time  over  an  extended  period. 
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in  the  USNO  Machine  Shop — Gary  Wieder,  Dave  Smith,  Tie 
Siemers,  and  John  Evans — for  fabricating  all  of  the  custom 
optomechanical  elements  of  dFTS2,  as  well  as  the  thermal 
enclosure.  We  also  thank  the  USNO  Astrometry  Department 
for  travel  support  and  salary  support  during  the  initial  stages  of 
this  observing  program,  and  thanks  also  go  to  the  USNO  Time 
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our  instrument  to  Kitt  Peak. 

This  research  has  made  use  of  the  SIMBAD  database, 
operated  at  CDS,  Strasbourg,  France;  NASA’s  Astrophysics 
Data  System;  and  the  SB9  catalog  of  Pourbaix  et  al.  (2004). 
Richard  O.  Gray  is  to  be  commended  for  making  his  SPEC¬ 
TRUM  codes  so  easy  to  install  and  use.  Thanks  also  go  to 
Farnoud  Kazemzadeh  for  a  critical  reading  of  the  manuscript, 
and  to  Frank  Fekel  and  David  Ramm  for  useful  email  discus¬ 
sions  regarding  the  finer  points  of  Spectroscopic  binary  analysis. 

Note  added  in  proof.  After  this  paper  went  to  press,  we 
were  alerted  to  the  work  of  Southworth  et  al.  2007  (A&A, 
467,  1215),  who  used  high-precision  photometry  to  derive  an 
inclination  angle  i  =  76:80  ±  O’:  1 0  for  the  beta  Aurigae  system. 
Using  this  value,  we  compute  alternative  results  for  Table  7: 
Mi  =  2.3644  ±  0.0045  M0  and  M2  =  2.3035  ±  0.0044  M0 
(formal  RV  errors)  or  M\  =  2.3644  ±  0.0055  M0  and  M2  = 
2.3035  ±  0.0055  M0  (scaled  RV  errors).  The  relative  mass 


errors  are  thus  0.19%  (formal)  or  0.24%  (scaled).  We  also 
overlooked  the  recent  publication  of  Konacki  et  al.  2010  (ApJ, 
719,  1293),  who  used  a  specialized  iodine  cell  technique  to 
achieve  excellent  radial  velocity  precision  on  l  Pegasi,  o>  Dra  and 
12  Bootis,  yielding  stellar  mass  values  with  relative  error  bars 
of  0.065%,  11%,  and  0.2%  respectively.  We  thank  the  authors 
of  both  papers  who  brought  these  references  to  our  attention. 
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